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ABSTRACT 



Using an SPH simulation of a star-forming region in a molecular cloud, we show 
that the emergence of a clump mass function (CMF) resembling the stellar initial mass 
function (IMF) is a ubiquitous feature of molecular cloud structure, but caution against 
its over-interpretation. We employ three different techniques to extract the clumps in 
this study. In the first two, we interpolate the SPH particle data to 2 and 3 dimensional 
grids before performing the clump-find, using position-position (PP) and position- 
position- velocity (PPV) information respectively. In the last technique, the clump- 
finding is performed on the SPH data directly, making use of the full 3 dimensional 
position information. Although the CMF is typically similar to that observed in regions 
of nearby star formation, the individual clumps and their masses are found to be 
unreliable since they depend strongly on the parameters and the method of the clump- 
finding. In particular we find that the resolution and orientation of the data make a 
significant difference to the resulting properties of the identified clumps in the PP and 
PPV cases. We conclude that making comparisons between a CMF and the stellar 
IMF should be done with caution, since the definition of a clump boundary, and hence 
the number of clumps and their properties, are arbitrary in the extraction method. 
This is especially true if molecular clouds are truly scale free. 
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Recent observations ( Motte et al, 

2006; N utter &; Ward-Thompson|(2007| |Alves et al 
Ward- Thompson et al.| 2007[ ) of dense gas cores in 



molecular clouds have shown their distribution resembles 
the stellar initial mass function (IMF). It is hoped that by 
studying these dense cores the earliest stages of star forma- 
tion can be probed and the origins of the IMF revealed. 

Motte et al.| ( |1998 ) have shown that the clump mass 
function (CMF) in p Oph can be described by a simi- 
lar power-law fit as the IMF. Above ^ 0.5 M© they find 
dn/dm oc is well- fit ted by an a = 2.5, while at lower 

masses they see a turn-over that can be fitted with a = 1.5. 
These values for a are broadly consistent with the usual fits 
to the IMF ( Kroupa]|2002 ). Similar results have been con- 
firmed by a number of authors for a variety of nearby star- 
forming regions, although the range of core masses found 
and the break mass of the power law fit can vary 



[stone et al1|2000l [2006| [Nutter Ward-Thompsonf|2007| 



John- 



Alves et al.[[2007[ [Testi & Sargent[[1998|. A more detailed 



account of the clump/core mass distribution can be found 



Ward-Thompson et al. (2007) 



One interpretation of the similarity between the CMF 
and IMF is that the masses of stars are controlled by the 
fragmentation of the cloud. [ Myers | ([2008| propose that for 
cores with a high density contrast to their surroundings, 
the resulting protostellar mass will be the gas mass whose 
free-fall time equals the core dispersal time due to outflows. 
Additionally [Swift Williams (2008) take a slightly differ- 



ent approach and argue that the resulting shape of an IMF 
formed from a Saltpeter-like CMF is robust against different 
core evolution scenarios. 

In this paper, we show that the shape of clump mass 
function seems to be a ubiquitous feature of molecular 
clouds, but the identified clumps and their masses are sen- 
sitive to the method used to find them. 

Molecular cloud (MC) structure as a whole is com- 
plex, filamentary and extends through all observable scales 
2000| ). It can also be thought of as hier- 
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archical: small dense cores are part of small clumps of gas 
which are themselves part of larger gas clumps. This con- 
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tinues over the full spatial range of observations (Bergin 



& Tafalla 2007). MCs also exhibit supersonic turbulence. 



Larson (1981) showed that the internal velocity dispersion 
<j{v) of a cloud is related to its size L by the relation 
(t{v) oc L~°"^^. Such a scale- free relation is a natural out- 



come of a turbulent cascade, with Brunt & Mac Low (2004) 



proposing that compressible shock dominated Burgers tur- 
bulence is most consistent with observations. 

As supersonic flows compress the gas in the molecu- 
lar cloud they create a complex filamentary density field, 
suggesting that the observed structure in MCs is a direct 
consequence of turbulent motions ( Klessen k, Burkert|2000 



[Kle ssen 2001). At the smallest scales of the hierarchy, the 
transonic dense cores formed by the turbulent cascade may 



be the precursors of star formation (Goodman et al.||1998 



Ward- Thompson et aL]|2007t . Indeed Pado an Nordlund 
(|2002|) have predicted that mass distribution of Jeans un- 
stable clumps created by super-Alvenic turbulence will re- 
semble the IMF, a prediction that is supported by the ob- 
servations of the CMFs discussed above. 

However there are still questions as to whether the CMF 
observations can unambiguously support the cloud frag- 
mentation picture. Typically, the observations must rely on 
two-dimensional column density maps (PP), which are lim- 
ited by resolution and sensitivity, and the observed features 
in the column density may be contaminated by integrated 
emission along the line-of-sight. Three-dimensional informa- 
tion is only available by obtaining position-position- velocity 
(PPV) cubes of the emission from a particular molecular 
tracer. Converting between emission and column density re- 
lies on many assumptions (Pineda et al. 2008) and account- 
ing for the properties of dust is itself challenging ( Schnee] 
et al.|[2006 ). Further, the observed emission is itself depen- 



dent on the entangled effects of density, temperature, molec- 
ular tracer properties, and the observational technique (Di 



Francesco et al.|2007| ). 



In this paper, we investigate how the differences in 
the form of the data - whether it be PP, PPV or full 
three-dimensional PPP - affects the resulting CMF. We ex- 
plore how the resolution the data and orientation of the 
cloud affect the mass functions profile. Since we are in- 
terested in comparing observational style clump-finds to a 
full three dimensional clump-find, the data must necessar- 
ily be synthetic. As such we use the results from a high- 
resolution smoothed particle hydrodynamics (SPH) simula- 
tion of driven turbulence. A full description of this simu- 
lation can be found in Section |2] To produce the PP and 
PPV style clump-finds the SPH data is interpolated to a 
grid, while for the PPP data the clump-find is done directly 
on the SPH data. The method used to extract the clumps is 



similar to that described by Williams et al. ( 1994), although 



there are some differences that can be found in Section |3] 



2 THE SIMULATION 

We use a three dimensional smooth particle hydrodynamics 
(SPH) code ( Monaghan 1 992), to simulate star formation in 
a collapsing gas cloud. The SPH code has variable smooth- 



ing lengths and artificial viscosity (Gingold & Monaghan 



1983 ) . Self- gravity is calculated via a binary tree ( Benz 1990 ) 
and sink particles ( Bate et 'aL]|1995 ) are used to model the 



Table 1. All the clump- finding presented in this paper were 
performed on a single SPH simulation. The initial conditions 
of our simulation are given here. The mass resolution is the 
minimum mass gravitational forces can be resolved for and is 
calculated via Mres ~ ^^^^ total Impart- The Jeans mass, free- 
fall time and sound speed are calculated from the average den- 
sity of the simulation before collapse using the formulae Mj = 
(47rp/3)-i/2(5iEcr/(2G/i))3/2, tff = (37r/(32Gp))i/2. 



Size 

Mass 

Particles 

Mass resolution 

Jeans mass 

Free-fall time 

Sound speed 



2 pc cube 
118.2 M0 
5000211 
2.36 X 10-3 M0 
5.7 M0 
2.1 X 10^ yrs 
0.18 kms-i 



star formation. Table [T] shows the initial conditions of the 
simulation. The molecular cloud has a mass of 118.2 Mq, 
roughly comparable to the masses of the 'cores' in p Oph 
(Mott e et al.|19 98). The cloud is modelled with 5 million 
SPH particles, which gives a mass resolution of 2.36 x 10~^ 
Mq ( Bate &: Burkert||1997 ). As such, we are able to resolve 
clumps with masses in the brown dwarf regime. The simula- 
tion was run on the SUPA Altix computer at The University 
of St Andrews. 

We ensure the velocity and density fields of the simula- 
tion are self-consistent as follows. A turbulent velocity grid, 
consistent with a Larson velocity dispersion of cr oc r°"^, is 



generated according to |Dubinski et al. (1995) and Myers & 
Gammie| (|1999). The SPH particle velocities are interpolated 



from this grid. The grid is continually re-imposed on the par- 
ticles as they are allowed to relax without self-gravity. The 
driving pattern is not fixed, but rather evolves between two 
realisations on a crossing time. After half a free fall time the 
probability distribution function has a lognormal shape as is 



expected for an isothermal gas without self-gravity, ( Klessen 
&: Burkert|2000 ), and the turbulence has power on all scales 



The turbulence has a r.m.s velocity of Mach 4.25 which cor- 
responds to the simulated gas cloud being just kinetically 
supported (Efc = |Ep|). 

When self- gravity is switched on, we stop the turbulent 
driving, allowing the supersonic motions to decay. |Mac Low| 
|et al. (1998) have shown that turbulence typically decays 
on a crossing time, however, the turbulence in our simula- 
tion remains supersonic throughout due to the conversion 
of potential energy to kinetic energy as the cloud contracts. 
Star formation, as modelled by the formation of the first 
sink particle, occurs after one free-fall time 2.1 x 10^ yrs. At 
this point the molecular cloud has the filamentary structure 
typically seen in observations (e.g. Bergin Tafalla||2007 ). 
We do not follow the evolution of the cloud beyond the for- 
mation of the first sink (that is, the runaway collapse of the 
first core). Rather, we use a snapshot of the cloud structure 
immediately before the formation of the first sink, as the in- 
put data for our clump-finding comparison. Figure [l] shows 
a column density projection of the cloud at this point in the 
evolution. 

A polytropic equation of state is used to mimic the evo- 
lution between the optically thin and optically thick regimes. 
Below densities of 10~^^ gcm""^ the value of the exponent 
gamma is initially set to 7 = 1 to mimic dust cooling. At 
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Figure 1. The column density of the data file used for the 
clumpfinding comparison. The density scale is logarithmic and 
runs from 0.005 gcm"-^ (black) to 1 gcm~^ (white). The region 
shown is 1 pc by 1 pc large and contains 79.5 Mq. 



gas densities of 10""^^ gcm~^ it changes to 7 = 7/5 where 
the gas becomes optically thick. Finally at densities greater 
than 10~^ gcm~^ it is set back to 7 = 1 to model the dis- 
sociation of H2 . It is during this last isothermal regime that 
sinks are allowed to form. 



3 CLUMPFINDING METHODS 

First, we wish to stress the nomenclature used in this paper. 



Williams et al. (2000) were among the first to try make a 



clear definition for the terms 'clump' and 'core' and Bergin 
|fc Tafalla (2007) have refined this. They define a clump as 
having a size 0.3 — 3 pc and containing 50 — 500 M0 and 
define a core as being around 0.03 — 0.2 pc and having a mass 
of 0.5 — 5 M0. For this paper we shall simply use the term 
clump to refer to any object that is over dense compared 
to its surroundings and has been identified by our clump- 
finding algorithm (that is, appears when the cloud is viewed 
in the chosen contours). By the above definition however, 
most of our objects are technically cores. 

The basic approach of the clump-finding algorithm is 
using a 'friend-of- friends' algorithm start from the point 
with the highest density, and move to lower densities, as- 
signing the cells/particles to the same clumps if they are 
neighbours. The first position will be assigned as the head 
cell of a clump, if the next highest emitting region is a neigh- 
bour of the previous one then it is assigned to the same 
clump, if not it is the head cell of a new clump. Contour 
levels are used to determine whether a clump is a significant 
feature or merely a bump on a larger clump. A similar tech- 



nique has been used by other authors (e.g. Klessen & Burk- 
ert|2000 ). If the heads of two neighbouring clumps and their 



but its head is on a lower density contour then the smaller 
clump is merged with the larger one. This process contin- 
ues until all the points have been considered. Although this 
method of clump-finding is similar to the publicly available 
CLUMPFIND algorithm (iWilliams et aL||1994|), we do in- 



clude a modification: the clumps are defined down to either 
a minimum density, or the lowest contour which it shares 
with a neighbouring clump. As such the mass in the clumps 
is typically much less than the mass in the cloud. By doing 
this, we are able to mimic the background subtraction that 
is normally performed on observational data sets. The def- 
inition of a 'clump' here is also more in keeping with those 
studied by |Motte et al.| ([1998 ) and P adoan et al.| ( |2007| . 

As Pineda et al (in preparation) have shown, clump- 
finding routines are very sensitive to the positions and num- 
ber of contour levels. In this study we ensured that the con- 
tour levels were static in all the cases considered here. The 
lower density sensitivity was usually set at a column density 
of 0.02 gcm~^, we investigate the effect of varying minimum 
density in Section [5] 

In this study, we apply our clump finding algorithm 
to three different styles of data set: position-position (PP), 
position-position-velocity (PPV) and position-position- 
position (PPP). The first two (PP and PPV) represent the 
type of data that one can obtain from molecular cloud sur- 
veys. In these cases the 'density' on which the clump-find 
is performed is the column density projected along the line 
of sight. In the PPV case, one simply has the projected col- 
umn density per velocity channel. We obtain these maps by 
interpolating the SPH particles to a grid of given resolution. 
Rather than interpolating to a full PPP grid and then inte- 
grating along the required line of sight, it is more efficient 
to integrate the spherically symmetric SPH smoothing ker- 
nel along the line of sight, and use the resulting 2D kernel 
to perform the interpolation. The full PPP clump-find is 
performed directly on the SPH particle data, although the 
process of clump-finding is exactly the same. 



4 TWO-DIMENSIONAL (PP) 
CLUMPFINDING 

As mentioned in Section [l] there have been a number of 
attempts to connect the structure of the gas in dense star- 
forming regions to the IMF. The majority of the work to date 
has produced results based on 2D maps of these regions, ei- 
ther using dust continuum emission 



Johnstone et al. 



2006) or extinction mapping (e.g|Lada et al.|| 200j| to ob 



intersection are on the same density contour the clumps are 
merged. If a clump shares neighbouring cells with another 



tain the column density distribution in the plane of the sky. 
As such, we start our clump-finding analysis on position- 
position column density maps of our simulated molecular 
cloud. 

Figure [2] shows the PP clumpfind we adopt as the stan- 
dard which we make comparisons against. The density is 
shown as a column density projection of the clumps after 
they have been background subtracted, the same scale is 
used as in Figure [l] for ease of comparison and to allow the 
clumps to be seen clearly, though no clump has a density 
lower than the minimum value of 0.02 gcm~^. In the regions 
where no clump has been assigned the lowest colour on the 
density scale is used. Table |2] shows the parameters used to 
obtain the PP standard clumpfind. 
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Table 3. The recovered clumps at different resolutions for clump- 
find on PP data. The * denotes our fiducial case and the resolution 
is given in grid cells. 



Resolution 


Spatial Resolution 


No. Clumps 


Total Mass 


50 X 50 


4, 125 AU 


12 


26.04 M0 


100 X 100 


2, 063 AU 


17 


18.62 M© 


200 X 200 * 


1,031 AU 


28 


16.35 M0 


400 X 400 


516 AU 


31 


15.71 M© 



-0.4 



-0.2 



0.2 



0.4 



Figure 2. The standard PP clump-find projected in the xy plane, 
crosses denote the centre of the clumps. This clump-find used 
data from a 200 x 200 grid. The scale is shown in parsecs and 
the colours represent column densities in the range 0.02 gcm~^ 
(blue) to 1 gcm~^ (yellow) at logarithmic intervals. 



Table 2. The parameters used to find our standard PP clumpfind 
shown in Figure [2] The resolutions grid cells correspond to a 
physical length of 1, 031 AU. 



Minimum Density 
Projection 
Resolution 
Spatial Resolution 
Density Countours 



0.02 gcm-'^ 
xy plane 
200 X 200 grid cells 
1,031 AU 
5 per magnitude 



4.1 Resolution 

Our analysis of the PP clump-finding procedure starts with 
examining a fundamental difficulty in interpreting the ob- 
served clump properties: would the results differ if the ob- 
served cloud was further away? We address this problem by 
changing the resolution of the grid onto which we interpo- 
late our SPH data. Perhaps unsurprisingly, we find that the 
number of clumps identified in the data varies with the res- 
olution of the grid. When the same physical region of the 
SPH data is interpolated to a grid with smaller elements, 
the clump-finding procedure yields more clumps. Figure [3] 
shows the clumps extracted by the algorithm from 4 data 
grids of the following sizes, 50 x 50, 100 x 100, 200 x 200 and 
400 x 400 grid cells . Table [3] shows the number of clumps 
found in each case and the spatial resolutions in physical 
units. Note that these resolutions are comparable to current 
observational surveys of nearby star-forming regions. Figure 
^ shows the clump mass functions. 

As the resolution of the clump-finding input increases 
the slope of the high mass end steepens due to more 
intermediate-mass objects being found and the biggest ob- 
jects being broken into smaller units. Additionally the break 
position between the two regimes decreases as the resolution 




Figure 4. The variation of the cumulative number clump mass 
function with resolution. The clump-finds shown are from the 
dashed line 50 x 50 grid with a cell width of 4, 125 AU, dotted 
line 100 x 100 grid with a cell width of 2, 063 AU, the solid line 
200 X 200 grid with a ceh width of 1, 031 AU and the long dashed 
line 400 x 400 grid with a ceh width of 516 AU. 



increases. The slope of the low mass end does not appear to 
change significantly, however the limited number of clumps 
in this regime, particularly in the low resolution cases, make 
it hard to draw definite conclusions. 

The amount of mass found in clumps decreases as the 
number of clumps and the resolution increases. There are 
two reasons for this. First, the higher resolution cases al- 
low one to resolve peaks in the structure that are blended 
at lower resolution. Secondly, we define our clumps down 
to the lowest contour that is shared with a neighbour. As 
such, two newly resolved smaller clumps that were previ- 
ously one larger clump, will have less mass than their larger 
counterpart. Note that in the standard CLUMPFIND algo- 
rithm, the total mass would be the same, since the clumps 
fill the entire surface/ volume of the data grid. However a 
similar analysis performed with CLUMPFIND on observa- 
tional data may still yield a similar result to what we find 
here, since the background is usually subtracted before the 
clump-finding is applied. 
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Figure 3. Four clumpfinds carried out on the central region where star formation wih occur with varying resolution levels. Top left A 
50 X 50 grid with a cell width of 4, 125 AU, top right a 100 x 100 grid with a cell width of 2, 063 AU, bottom left a 200 x 200 grid with a 
cell width of 1, 031 AU and bottom right a 400 x 400 grid with a cell width of 516 AU. The scale is shown in parsecs and the grayscale 
represent column densities in the range 0.02 gcm~^ (grey) to 1 gcm~^ (black) at logarithmic intervals. Crosses show the centre of the 
clumps. 



4.2 Orientation 

A further difficulty in interpreting the results of PP clump- 
finding is that a three dimensional structure has been ar- 
bitrarily projected into two dimensions, hence, it is unclear 
whether a clump is one object, or a superposition of sev- 



eral objects along the line of sight ( Ballesteros-Paredes 
Mac Low 2002 ) . Figure [s] shows position along the line of 
sight and the real 3D density of the material contributing to 
a typical dense clump in the standard PP clump-find. The 
material is too widely distributed to be a single core. More- 
over, some of this gas lies at densities (3D) which actually 
sufficiently preclude it belonging to any dense core. 

Figure [6] shows that line of sight blending causes differ- 
ent clumps to be found if the three dimensional data set is 



Table 4. The properties of the clumps found in different orien- 
tations. The starred case corresponds to our standard PP clump- 
find. 



Projection 


No. Clumps 


Total Mass (M©) 


xy * 


28 


16.35 


xz 


10 


27.33 



projected along a different axes. Table [3] shows that different 
orientations result in different clump properties. This is in 



agreement with previous work by Ostriker et al. (2001 ) 
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Figure 5. The intrinsic 3D density profile along the z axes of the 
material belonging to a typical clump in the standard PP clump- 
find. About half the material assigned to the clump is truly at 
high 3D densities but the rest is contamination from material 
along the line of sight. 




Figure 7. The cumulative number clump mass functions found 
from PP data with three different lower column density limits; 
solid line 0.01 gcm~^, dashed line 0.05 gcm~^ and dotted line 
0.1 gcm~^. 




Figure 6. A clump-find with identical parameters to our stan- 
dard PP clump-find carried out on the simulated data projected 
into the xz plane instead of the xy plane. The scale is shown in 
parsecs and the grayscale represent column densities in the range 
0.02 gcm~^ (grey) to 1 gcm~^ (black) at logarithmic intervals. 
Crosses show the centre of the clumps 



4.3 Density Cut 

The result of a clump-find is also dependent on the den- 
sity range of the data-set. Figure [7| shows clump mass func- 
tions obtained from a basic PP data-set with three different 
lower-density sensitivities. As the density is lowered more 



clumps are discovered and a greater fraction of the struc- 
ture is probed. 

With a restricted density range only the high mass end 
of the CMF is found. Then when the sensitivity is lowered 
further the CMF resembles two power laws. 

If the PP data spans a narrower density range the re- 
sulting clumps are more distinct. Due to their decreased size, 
they contain less mass, in other words only the dense cen- 
ters of the clumps are now included. The second panel of 
Figure [8] shows what the standard PP clump- find becomes 
if applied to gas with a number density above 10^ cm^. 



5 POSITION-POSITION- VELOCITY (PPV) 
CLUMPFINDING 

The second type of data set we consider in our analysis is 
position-position column density maps with the radial ve- 
locity used as an additional dimension. This data set should 
suffer less from line of sight blending as it is likely that dif- 
ferent clumps will be travelling at different velocities. The 
width of the velocity channels used determine the velocity 
resolution of the clumpfind, in our analysis we assign 16 
velocity bins along the z-axes of our data cube which corre- 
sponds to a resolution of 0.33 kms~^. Due to the increased 
number of grid cells a PPV clumpfind uses, it has an intrin- 
sically higher resolution than a PP clumpfind. Panel three 
of Figure |8] shows the PPV clumpfind, and Table [5] compares 
its properties to the standard PP case and to an additional 
PPV case with 40 velocity bins. 

In the PPV case, clump masses are reduced as mate- 
rial previously included is rejected due to its velocity. The 
positions and extent of the clumps also change relative to 
the PP clumpfind. Changing the size of the velocity bins 
also changes the mass and number of clumps found. In the 
case shown here, when the bin size is decreased and hence 
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Figure 8. Clumpfinds using various data-sets. Top left the standard PP clump-find. Top right applied to a data set above an intrinsic 
3D density of lO'^ cm~^. Bottom left using the PPV data-set with 16 velocity bins, corresponding to a velocity resolution of 0.33 kms"-*^. 
Bottom right the 3D clumpfind on the raw SPH data. The scale is shown in parsecs and the grayscale represent column densities in the 
range 0.02 gcm~^ (grey) to 1 gcm~^ (black) at logarithmic intervals. Crosses show the centre of the clumps. 



Velocity Bins 


Resolution 


No. Clumps 


Total Mass 


0* 




30 


16.37 M0 


16 


0.33 kms-i 


37 


8.90 Mo 


40 


0.13 kms-i 


33 


4.50 Mo 



Table 5. The properties of the clumps using the PPV data-set 
compared to the standard PP clump-find. The 16 velocity bins 
have a width of 0.33 kms"-*^ and the 40 velocity bins a width of 
0.13 kms-i 



the resolution is increased, there is less mass assigned to 
the clumps. A side effect of this process is that some of the 
clumps masses decrease to the point that they no longer 



meet the resolution requirement and are therefore not in- 
cluded in our analysis. 

As regards blending, if we consider the clump found 
using the PP method shown in Figure [5] Figure [5] shows 
that more of the mass is now contained in a smaller region. 
However, the density profile suggests that the clump found 
still contains spurious material that happens to be mov- 
ing towards the centre of the region with a similar velocity. 
Nonetheless, the contribution from material at low density 
along the line of sight has been reduced. 

It should be noted that the analysis presented here for 
the PPV method is a 'best case scenario' as no molecu- 
lar tracer can actually probe the entire line of sight from 
the low density outer layers to the dense core centre. This 
means that the velocity information is likely to be incom- 
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z(pc) 

Figure 9. A projection of the real 3D density profile of the clump 
shown in the lower panel of Fig [s] when found using PPV data 
with 16 velocity bins. 



plete. Further, temperature gradients in the core centres will 
also affect the line profiles. An additional caveat of the PPV 
method is that if the clump is undergoing collapse or ex- 
pansion then it could be split into two objects depending on 
the velocity resolution. However this is still a useful observa- 
tional technique, particularly when using multiple molecular 
tracers (e.g. |Andre et al.||2007 ). 



6 THREE-DIMENSIONAL (PPP/SPH) 
CLUMPFINDING 

The final data set we consider is that of our full simulated 
three-dimensional structure. While this is impossible from 
observations, it is still a useful comparison as the true struc- 
ture of star forming regions is, of course, three-dimensional. 
The same basic method is used to find the clumps as before. 
This is possible as in SPH forces are calculated using a hi- 
erarchical tree, consequently a neighbours list exists which 
allows the densities of neighbouring particles to be compared 
in a similar manner to grid cells. 

However, unlike the PP and PPV clump-finds the reso- 
lution is no longer simply the size a grid. In SPH simulations 
the mass resolution varies spatially, high density regions are 
better resolved than low ones, therefore our three dimen- 
sional clumpfind also shares this property. This effectively 
results in our SPH run having the highest resolution of all 
the clumpfinding methods. 

The bottom right panel of Fig. [S] shows the result of 
the clump-find on the 3D SPH data. The number of clumps 
has increased and the total mass decreased relative to our 
previous results. There are two reasons for this. First, due to 
our improved resolution the number of clumps increases as 
shown in Section [XT] Secondly the use of 3D data eliminates 
the blending along the line of sight previously encountered in 
Section [4.2| so more clumps are found with smaller masses. 
Effectively we are now detecting a finer level of structure 



Method 


No. Clumps 


Total Mass (M©) 


PP 


28 


16.35 


PP high density 


19 


10.51 


PPV 16 bins 


37 


8.90 


sph 


72 


4.50 



Table 6. A comparison of the clump-find results. 

than in the previous data-sets, namely the very dense core 
centers. 



7 DISCUSSION OF GLOBAL PROPERTIES 

7.1 The Clumps 

Table |6] compares the properties of the clump-finds. Al- 
though the general outline of the dense region where the 
clumps are identified remains constant, the number, masses 
and assigned boundaries of the clumps varies with extrac- 
tion method. Generally the number of clumps increases and 
the total mass found in clumps decreases with the precision 
of the method. 

Applying a density cut to a data-set produces more dis- 
tinct clumps, but is of course completely arbitrary. The PPV 
clump distribution is a better match to the high-resolution 
3D clumpfind than the PP clumps. The peak positions iden- 
tified as the heads of clumps from the PPV data corre- 
lates well with the highest density structures seen in the 
sph case. The overall mass distribution and the clump prop- 
erties, however, differ. Each data-set appears to highlight a 
different scale of the structure in our simulation. 

Defining the clump boundaries is further complicated 
by arbitrary decisions that have to be made in the 
clumpfinding process. Is there a minimum density required 
for a clump? How many density contour levels need to be 
set? Should some minimum contrast to the background be 
required? All these issues add to the difficulties presented 
here in comparing clumps found from different data-sets. 
These challenges should be kept in mind when comparing 
the properties of cores recorded in the literature, particu- 
larly when contrasting masses and sizes. 

7.2 The Clump Mass Function 

Despite the differences in the clump properties, the clump 
mass functions show marked similarities. Figure [To] shows 
the cumulative clump mass functions derived from the PP, 
PP with density cut, PPV and PPP (SPH) clump-finds and 
Table [7| shows the gradients of the power laws which best fit 
them. 

As the clump properties differ each of the CMF's has a 
different magnitude, size range and break point but the over- 
all shape is always broadly consistent with the stellar initial 
mass function (| Kroupa|2002 | . When plotted logarithmically 
this corresponds to a Saltpeter slope of F = —1.3 with a 
break to second power law of F = —0.3 at lower masses. 

At the low mass end, with the exception of the PP high 
density case, all the gradients are within 0.15 of the Kroupa 
value and within 0.08 of each other. The PP case where 
the density range was restricted to high values resulted in 



The structure of molecular clouds and the universality of the clump mass function 9 




S 10 



I — ' — 




J I I I I i\n I I I 



0.01 



10 




J I I 1. 1 1 1 1 1 I I I I I I 



Figure 10. The cumulative number clump mass functions from the previous clumpfinds. Top left our standard PP clump-find, top right 
high density PP data, bottom left the PPV clump-find and bottom right the PPP (SPH) method. The dashed lines are power law fits 
whose gradients are shown in Table [T] 



Method 


Low Mass Gradient 


High Mass Gradient 


PP 


-0.37 


-1.15 


PP high density 


-0.16 


-1.35 


PPV 16 bins 


-0.39 


-1.10 


sph 


-0.45 


-1.99 



Table 7. The best fit power law gradients for the cumulative 
CMF's shown in Figure [lO] 



incomplete sampling of the low-mass end, which leads to a 
flattening of the low-mass gradient. At the high mass end the 
two-dimensional methods have gradients within 0.25 of each 
other, poisson noise makes these fits more uncertain but the 
agreement is still good. The best fit to the high mass end 
of the SPH cumulative clumps mass function, however, is 
steeper than the two dimensional cases. It is unclear whether 
this is a genuine consequence of the 3D data set or due to 
the increasing resolution of the SPH method with density 
preferentially identifying smaller denser structures. 

Further evidence that the CMF may have a universal 
shape can be seen in our earlier discussions of resolution 
(Figure [4| and density range (Figure |7|. Once again in these 
cases we see that although the properties of the clumps are 
changing, the profile of the CMF remains broadly consistent 
with the steflar IMF. 



These results imply that the shape of the CMF results 
naturally from breaking up the underlying structure of a 
molecular cloud in a consistent manner. Therefore, the emer- 
gence of a rational CMF from a data set is no guarantee 
that the identified clumps are unique objects. The structure 
could have been divided up using a differing, perfectly jus- 
tifiable, method which results in different clump positions 
and sizes but yields a similar CMF profile. This is reminis- 
cent of the observations discussed in the introduction, all 
have similar CMFs but the characteristic masses of clumps 
differ. Indeed, due to difficulties like those discussed here. 



Rosolowsky et al. (2008) has proposed that dendrograms 
might ultimately prove a better way to identify structure 
than the traditional clumps. Still, many of the issues raised 
here are likely to be relevant using dendrograms or any other 
clump finding algorithm. 

7.3 Universality of the CMF 

Despite variations in clump properties, the slope of the 
clump mass function is similar in all the cases investigated 
here, with the possible exception of the high mass end in 
the SPH clump mass function. The CMF retains this shape 
when extracted using different clump-finding methods and 
parameters. 



As many studies have shown (e.g. Vazquez- Semadeni 
eraL][T995l [Ballesteros-Paredes et aLl|1999| |Klessen||2001| ), 
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turbulence generates the characteristic multi-scalar struc- 



ture of molecular clouds. Elmegreen (2002) goes so far as to 
describe the structure of the interstellar medium as a con- 
tinuous fractal distribution. Observationally, Padoan et al.| 
(2003) show that the structure of the Taurus and Perseus 

3 pc. Given 



molecular clouds is hierarchical at scales of 0.3- 



the universality of turbulence (Heyer & Brunt 2004), it is 



perhaps unsurprising that the observed CMFs exhibit sim- 
ilarities between different star forming regions, and there 
exists in the literature a number of models for describing 
how these clump s can form ([Fleck Jr. [1982 ElmegreenI 19931 
Padoan 1995; Padoan, Nordlund & Jones|1997i Myers 2000; 



Klessen 



2001| [Padoan Nordlund||2002t . The similarities 



between the CMF and the IMF then naturally leads to the 
sugge stion that the two are related ([Motte, Andre Neri 



19981 |Lada Muench, Rathborne, Alves LombardipOOSj 
Myers|20 08t" 

However, our findings suggest that the mapping be- 
tween the CMF and IMF is not straightforward. Analysing 
fractal-like structure with any clump-finding algorithm will 
generate a clump mass function. The details of the clump 
properties, especially the mass, will depend on the type of 
algorithm and the parameters chosen. Further, the clumps 
are not uniquely defined. This has also been found obser- 
vationally. [Johnstone et al.l (|2000|) and [Motte et al.[ ([1998) 



both surveyed the p Ophiuchi and derived clump mass func- 
tions which were in agreement with each other. However, as 
Johnstone et al. (2000) comment, the masses, numbers and 



positions of the clumps found do not correspond between 
the surveys. Due to the variance of clump properties with 
identification method, it is unclear how best to link star-less 
cores to protostars in a quantitative manner. 



8 CONCLUSIONS 

We have shown that when a CLUMPFIND style algorithm 
is applied to different types of data-sets from the same un- 
derlying simulation, clump characteristics vary. This is due 
to a number of factors. Details of the observational maps 
or the simulation obviously play their part, particularly pa- 
rameters such as resolution, orientation and density range. 
The more levels of structure resolved within a data-set the 
greater the number of clumps that can be found. Blending 
of structure along the line of sight decreases the number of 
clumps found in 2D clumpfinds and increases clumps masses. 
A significant amount of material which is in reality not con- 
nected to a high density core can be artificially included. 
Applying a high density cut to PP data yields more distinct 
objects but can bias the resulting mass function. 

Using velocity as a third dimension improves accuracy 
and provides a closer match to the true 3D structure than 
just using PP data. The clump boundaries, however, still 
differ and less mass is assigned to clumps. 

Qualitatively the profile of the clump mass function 
always resembles the Stellar IMF, regardless of the data- 
set used. However, the quantitative values of the CMF are 
dependent on the extraction method, and therefore causal 
comparisons to the IMF should be made with caution. 
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